Symmetry-Breaking and Symmetry-Restoring Dynamics of a Mixture of Bose-Einstein 

Condensates in a Double Well 
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We study the coherent nonlinear tunneling dynamics of a binary mixture of Bose-Einstein condensates in a 
double-well potential. We demonstrate the existence of a new type of mode associated with the "swapping" of 
the two species in the two wells of the potential. In contrast to the symmetry breaking macroscopic quantum 
self-trapping (MQST) solutions, the swapping modes correspond to the tunneling dynamics that preserves the 
symmetry of the double well potential. As a consequence of two distinct types of broken symmetry MQST 
phases where the two species localize in the different potential welils or coexist in the same well, the corre- 
sponding symmetry restoring swapping modes result in dynamics where the the two species either avoid or 
chase each other. In view of the possibility to control the interaction between the species, the binary mixture 
offers a very robust system to observe these novel effects as well as the phenomena of Josephson oscillations 
and pi-modes. 



I. INTRODUCTION 



Ultracold laboratories have had great success in creating 



Bose-Einstein condensates (BECs) in a variety of atomic 
gases such as Rubidium (Rb), Lithium (Li), Sodium (Na) and 
Ytterbium (Yb). These quantum fluids exist in various iso- 
topic forms as well as in different hyperfine states. The rapid 
pace of development in this field has led to condensates which 
are robust and relatively easy to manipulate experimentally. 
In particular, the tunability of inter-species and intra-species 
interactions 101 via magnetic and optical Feshbach resonances 
makes the BEC mixture a very attractive candidate for explor- 
ing new phenomena involving quantum coherence and non- 
linearity in a multicomponent system. 

The subject of this paper is to investigate the tunneling dy- 
namics of a binary mixture of BECs in a double well poten- 
tial. A single species of BEC in a double well is called a 
bosonic Josephson junction (BJJ), since it is a bosonic analog 
of the well known superconducting Josephson junction. In ad- 
dition to Josephson oscillations (JO), the BJJ exhibits various 
novel phenomena such as 7r-modes and macroscopic quantum 
self-trapping (MQST), as predicted theoretically |j3i|41]. In the 
JO and the 7r-modes, the condensate oscillates symmetrically 
about the two wells of the potential. In contrast to this, the 
MQST dynamics represents a broken symmetry phase as the 
tunneling solutions exhibit population imbalance between the 
two wells of the potential. These various features have been 
observed experimentally Isj]. Our motivation is to explore 
whether new phenomena arise when there are two interacting 
condensates trapped in a symmetric double well. 

Although our formulation and results are valid for a vari- 
ety of BEC mixtures, our main focus here is the Rb family 
of two isotopes, namely the mixture of *^Rb and ^^Rb, mo- 
tivated by the experimental setup at JILA|0]. The scattering 
length of ^^Rb is known to be 100 atomic units while the inter- 
species scattering length is 213 atomic units. In experiments, 
the scattering length of *^Rb can be tuned using the Feshbach 
resonance method 

The ability to tune the scattering length of one of the species 
makes this mixture of isotopes an ideal candidate for studying 



the coupled BJJ system. First, it opens up the possibility of 
exploring the parameter space where the Rb 85-85 scatter- 
ing length is equal to the Rb 87-87 scattering length. As will 
be discussed below, this symmetric parameter regime simpli- 
fies the theoretical analysis of the system and also captures 
most of the new phenomena that underlie the dynamics of the 
binary mixture. Furthermore, the tunability of the ®^Rb scat- 
tering length can be exploited to study a unique possibility 
where one of the species has a negative scattering length, a 
case which strongly favors the 7r-mode oscillations that have 
not been observed so far. 

In our exploration of nonlinear tunneling dynamics of cou- 
pled BJJ systems, the MQST states are found to be of two 
types. In the broken-symmetry MQST state, the two com- 
ponents may localize in different wells resulting in a phase 
separation or they may localize in the same well and hence 
coexist. By varying the parameters such as initial conditions, 
the phase-separated broken-symmetry MQST states can be 
transformed to a symmetry-restoring phase where the species 
continually "avoid" each other by swapping places between 
the two wells. In other words, if the dynamics is initiated with 
both species in the same potential well, the sustained tunnel- 
ing oscillations are seen where the two species swap places 
between the well one and the well two. From the coexisting 
MQST phase, one can achieve symmetry restoring swapping 
dynamics by initiating the dynamics with two species in the 
separate wells. In this case, the emergence of the swapping 
modes can be interpreted as a phase where the two species 
"chase" each other 

The paper is organized as follows. In section Ull we discuss 
the model and use the two-mode approximation to the Gross- 
Pitaevskii (GP) equation to map it to a system of two coupled 
pendulums with momentum-dependent lengths and coupling. 
Section Hill discusses the stationary solutions and their stabil- 
ity. These results enable us to look for various qualitatively 
different effects without actually solving the GP equations. 
Section|IV]describes the numerical solutions of the GP equa- 
tions as various parameters of the system are tuned. Although 
we have explored the multi-dimensional parameter space, the 
novelties attributed to the binary mixture in a double well trap 
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are presented in a restricted parameter space where the scatter- 
ing lengths of the two species are equal. Additionally, in our 
numerical results described here, we fix the ratio of Rb 87-87 
interaction to Rb 85-87 interaction to be 2.13. This restricted 
parameter space is accessible in the JILA setup and provides 
a simple means to describe various highlights of the mixture 
dynamics. Section [Vl provides additional details of the JILA 
setup relevant for our investigation. A summary is given in 
section 



II. TWO-MODE GP EQUATION FOR THE BINARY 
MIXTURE 

In the semiclassical regime where the fluctuations around 
the mean values are small, the two-component BEC is de- 
scribed by the following coupled GP equations for the two 
condensate wave functions $/(a;,t), with / = a, 6 represent- 
ing the two species in the mixture. 



2ma 

Here, mi, Vi and gi = Anfi^ai/mi, denote respectively, the 
mass, the trapping potential and the intra-atomic interaction of 
each species, with a; as the corresponding scattering length. 
gab = 27rfi,^(l/ma + l/mb)aab is the inter-species interac- 
tion, where Gab is the corresponding scattering length. For the 
JILA experiment, in view of the tight confinement of the con- 
densate transverse to the trap, it is sufficient to consider the 
corresponding one-dimensional GPE equations. 

The condensate wave functions satisfy the normahzation 
conditions. 



Ni 



(1) 



The total number of atoms in the mixture is iV = Na + N},. 
In the weakly linked limit, the dynamical oscillations of the 
two-component BEC can be described by two wave functions 
representing the condensate in each trap labeled by fc = 1,2, 
with the spatial and the temporal contribution factored as fol- 
lows: 



(2) 



The localized spatial modes Xk"^ i^) are computed as sums 
and differences of the symmetric and antisymmetric solutions 
of the time-independent, coupled GP equations. 

To derive the equations of motion in the two-mode approx- 
imation, we introduce zi{t), the population imbalance, and 
4>i{t), the relative phase of of species I between the left and 
right sides of the double well potential. 



where V'fc^(t) = \^k'it)\ exp iO^. are the time-dependent 
coefficients in the two-mode equations .Substituting equa- 
tions ^ into the coupled GP equations and integrating over 
the spatial degrees of freedom, we obtain the following four 
coupled, nonlinear, ordinary differential equations which we 
refer to as the "two-mode" model. 



it) = io[ 



(4) 



a(i) 



Za = -Ka\J\ - Z'l sin (/la 

Zb = -KbJl- Zl sin^fo 



= J^afaZa + KbfbZb + K 
h = ^bfbZb + AabfaZa + Kb 



Zg 

Zb 



(5) 
(6) 

cos 0a (7) 

COS (/)(,. (8) 



where Zi — zi/ fi. In the above equations, // ~ Ni/N de- 
notes the fraction of atoms of species /, while Ki and A/ are 
given by 



Ka - Xa - 2/a Ca^r^^ COS </.a + /bDaM/ 1-^2 COS 



Kb 
Aa 



Kb - 2fbCb y 1 - COS (j)b + faDba \/l ~ Zl COS 0a 

Aa +Ca 
Ab + Cf, 



In the above, the space and time-independent parameters Ki, 
A;, Aah, Ci and Dab can be expressed in terms of various mi- 
croscopic parameters that appear in GP equation and the lo- 
calized modes, xi (s^) and their overlap (integrated over spa- 
tial degrees of freedom). The explicit expressions for these 
parameters are given in the Appendix. 

The parameters Ki describe the tunneling amplitude while 
A; is related to the corresponding scattering length of the 
species. The parameters C; and Dab have their origin in the 
overlaps between the spatial modes xi and X2, and are ex- 
pected to be small in the weak tunneling limit. These overlaps 
modify the bare parameters denoted by the interaction A; and 
the tunneling /v;. Consequently, we have a variable tunneling 
model, since the tunneling parameters Ki depend explicitly 
on the dynamical variables Zi and 0/. 

In our analysis, we will mostly restrict ourselves to the 
case where the two species are equally populated, namely 
/a = /h = 1/2. In this case, the above system of equations 
(O can be viewed as the Hamilton equations in terms of the 
canonical variables Zi (momenta) and 0/ (co-ordinates), with 
the Hamiltonian given by the following form: 



n 



KaZl+KbZl- 



-2KabZaZb] 



l—a,b 



K,Jl 



Zf cos (j)i. 



(9) 



Zl it) = iM 



I |2 



\^ir)/N, 



(3) 



For the case where the overlap between the spatial modes 
Xi and X2 can be neglected, and the effective tunneling Ki 
can be replaced by its bare value Ki the above system can 
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be viewed as a coupled pair of non-rigid pendulums, with 
momentum-dependent lengths. The coupling between the 
pendulums is also momentum dependent. 

We parenthetically remark that this system can also be 
mapped to a pair of classical spins with Cartesian components 



Si 

si 



1 — Zf COS 0i 



1 - Zf sin 9i 



= Zi, 



so that (5')^ = 1. Thus the spin vector locates a point on 
the unit sphere given by polar angles 9i,<j>i, with Zi = cos 9i. 
The corresponding spin Hamiltonian, written in terms of bare 
variables, can be shown to be 



l—a,b 



The spin mapping provides an alternative means to visual- 
ize the effective interaction between the two species during 
the tunneling. If we ignore the spatial overlap integrals be- 
tween the localized modes in two wells, ( C; = 0, Dab = 
), the binary mixture of condensates in two-mode approxima- 
tion, maps to two Ising-like spins in a transverse magnetic 
field. The full two-mode variable tunneling feature induces 
XY-like spin interaction. 

In this paper, we find it convenient to exploit mapping to the 
coupled pendulums, for exploring tunneling dynamics in the 
coupled BJJ. Although we have explored the full two-mode 
variable tunneling model, we will only discuss the constant 
tunneling case ( Ki replaced by Ki and A; replaced by A;.), 
as the overlap integrals are small and the various novel effects 
of the mixture described here are found to be robust and unaf- 
fected by the variable tunnehng parameters. 



in. STATIONARY SOLUTIONS : FIXED POINTS 

The solutions of the coupled system are characterized by 
the interactions A;, the ratio of the tunneling amplitude for 
the two species, Ka/Kb which we denote by R as well as 
the initial phase difference — 0) and the initial popula- 
tion imbalance Zi{t = 0). In the multi- dimensional param- 
eter space the equilibrium or fixed-point solutions, in which 
the right-hand-sides of Eqs. (??)-(??) are zero, provide an ef- 
fective tool to classify different categories of behavior of the 
system. 

In general, these fixed-point equations are transcendental 
and have to be solved numerically. However, in the symmet- 
ric case where A^ = A;, = A, , Ka = Kb = K, the fixed 
point equations can be tackled analytically. Further, as can 
be seen from Eqs. (??)-(??), the parameter K can be elimi- 
nated in this case by rescaling t {t Kf) and redefining 
as Kx ~> Kx/ K. Our detailed analysis shows that this special 



case captures many relevant phenomena characterizing the bi- 
nary mixture in a double well. In this case, the fixed points be- 
long to two broad categories as stated below, resulting in two 
types of small amplitude oscillations about these two fixed 
points. It is important to note that this type of MQST phase 
does not exist in a BJJ with a single species. 

(I) Zero-mode Fixed Points (0^ 0* = 0) 



(1) zi = z*=o 

KD - -Z^ - ± |(A^^_A)| 

(II) TT-mode Fixed Points ( (; 



(1) K 



zt 







x/(Aab+A)"-4g^) 



(2) ZI - Z'l - ± (A^j,+A)| 

It should be noted that the mixed-mode Fixed Points, (</)* = 
and</.,*=7r),(Z: = Z*=0) 

are unstable for the restricted parameter regime we are 
considering here and hence will not be discussed. 

The small oscillations about the fixed point {Z* = 0, 0* = 
0) result in Zero-mode while small oscillations about {Z* = 
0, (/)* — tt) lead to 7r-mode. The oscillation frequencies are in 
the next subsection. 

The non-trivial fixed points {Z* 7^ ) result in solu- 
tions with population imbalance and lead to tunneling dynam- 
ics with macroscopic quantum self-trapping or the MQST. In 
view of the a — h symmetry, we have two sets of stationary so- 
lutions: Z* and — Z^), {x = a, 6) This suggests the possibility 
of modes where each species oscillates about the binary fixed 
points, going back and forth between the two wells. Unlike 
MQST, these modes will preserve the symmetry of the dou- 
ble well. However, in contrast to Zero-modes, these modes 
are non-linear and give rise to "swapping phase" that will be 
discussed later. 

The emergence of fixed points with opposite signs for the 
two species, (Z* = — in the Zero mode phase suggests 
that MQST in Zero-mode is accompanied by phase separation 
of the two species. In contrast, in the 7r-mode MQST phase 
, the two species could coexist in the same potential well as 
(Z* ~ Z^). Therefore, the fixed point equations suggest that 
TT-modes mimic attractive interaction between the two species. 

The onset from oscillatory to MQST phase corresponds to 
the values of the parameters where the non-trivial fixed points 
move from the complex to the real plane. Alternatively, the 
condition for the broken symmetry phase can be obtained by 
linear stability analysis of the fixed-point equations. This is 
discussed in the next sub-section. 

In the asymmetric case when A^ 7^ A^ the fixed points are 
obtained by solving the coupled transcendental equations: 



K 7* 

{-Vf "a 



1 



KhZl , 1 



{t.zi + KbZl) 

[Kzi+Kbz:) 



where p ~ 0(1) for 0* — 0(7r) and 0^ ~ 0(7r).. 





0, 
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Analogous to the symmetric case, both the Zero and the tt- 
mode solutions including those corresponding to MQST can 
be found numerically. As expected, for the MQST fixed points 



in the Zero-mode and Z* ^ in the 7r-mode 
and we do not have the permutation symmetry or the a — b 
symmetry. However, unlike the symmetric case, A'/s do not 
scale time t and the parameters and hence the ratio R 
emerges as a new parameter. 



Kg 

Kb 



A. Normal Modes: Linear Stability Analysis of Fixed Points 



Frequencies of small amplitude oscillations about about 
{Z* = 0, = 0) and (Z* = 0, = tt) respectively re- 
ferred to as the Zero-mode or the 7r-modes are given by 



where 



K 



faK 

fb^b, 




where p = for the Zero-mode , and p ^ 1 for the vr-mode. 
In the symmetric case, with Aa = At and fa = /b, the normal 
mode frequencies oja and oj^ simplify to. 



K{A±Aab)/2 
K{A±Aab)/2 



The condition for the instability of the fixed point is de- 
termined when one of the normal mode frequencies become 
complex. This gives rise to new fixed points where Z* ^ 
resulting in MQST phase where there is a population imbal- 
ance between the two wells of the double well potential for 
each species. The condition for the existence of MQST is 
given by. 



fafbA^b > 



A* A* 



(10) 



For the parameter values where both the Zero and the tt- 
modes coexist, 7r-mode frequencies are smaller than the Zero- 
mode frequencies. 

Figure 1 shows the values of Aa, Ab where the tunneling is 
governed by the Zero-mode and the 7r-mode. For A;, > 0, the 
regime where the 7r-modes exist is is small but finite. How- 
ever, by tuning A;, to negative values, the 7r-modes that have 
not been seen in earlier studies, can be observed. Variation 
with the parameter R, the tunneling ratio for the two species, 
leads to similar results, with the parameter space for the ex- 
istence of TT mode increasing slightly with R. The unshaded 
regime corresponds to MQST phase. 



FIG. 1: (color online) The upper(black;) and the lower(yellow) 
shaded regime corresponds to the parameter values for the existence 
of stable Zero-mode and 7r-mode with R = 1. 



numerical solutions also illustrate nonlinear modes, not de- 
scribed by the fixed point analysis. The fact that new features 
continue to exist in the nonlinear regime, assures their robust- 
ness. 

In our numerics, we set Aah = 2. 13 A and study the dy- 
namics for different values of A. These conditions can be 
achieved by first tuning the gj, via a Feshbach resonance so 
that Aa = Aft. The variation of A corresponds to varying 
the number of atoms in the double well trap. As already men- 
tioned, K can be eliminated by using t Kt and A ^ A/K. 
The dynamics is governed by A and the initial conditions: 
Za(0),Zfc(0), 0a(O) and (/.^(O). 

As we discuss below, tunneling solutions belong to three 
broad categories: 



(I) " Zero-phase Mode " , characterized by < (/); >= 



(II) "TT-phase Mode " characterized by < 0z >= tt 



(III) " Running-Phase Mode " characterized by < 0j > 
proportional to t 



IV. TUNNELING DYNAMICS WITH Aa = At 

We now describe numerical solution of the tunneling equa- 
tions. For small population imbalance, we confirm the dynam- 
ics predicted by the fixed points as discussed above. However, 



In the single species case, <(/);>= also corresponds 
Xo < Zi >= {]. However, as we discuss below, in a binary 
mixture, we can have < 0/ but < Z; >^ 0. This 
gives rise to a broken symmetry MQST phase in Zero-modes 
as well. 
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Zero-Modes 



that species move more sluggishly in 7r-mode compared to the 
Zero-mode as the Zero-mode frequencies are larger than those 
of the TT-mode. 



For A < « 1.77, and 0a (0) = 06 (0) = 0, and 
|Z((0)| << 1, both species execute small amplitude os- 
cillations (like oscillations of a non-rigid pendulum) with, 
< Zi{t) >= and < >= as shown in Fig. 2. Such 

modes exhibit quasiperiodic dynamics characterized by super- 
position of sinusoidal modes with two competing frequencies. 
As Zi (0) increases, we see large ampUtude non-sinusoidal os- 
cillations. Therefore, in spite of the repulsive interaction be- 
tween the two condensates, the two species execute a coherent 
oscillatory dynamics as expected from the zero-mode fixed 
point analysis described earlier. 





-0.2 
( A=.6, n- modes ) 



FIG. 2: (color online) Time series for Zi(lop), <j)i (middle) and phase 
portrait for A = 0.6 with initial conditions shown in the figure. The 
red and blue corresponds to the a and b species. 



B. TT-modes 
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FIG. 3: (color online) Same parameters as Fig. 2, the only exception 
being that 0j = 0) — n here. 



Symmetry Breaking and Phase Separation: MQST in 
Zero-Mode 



Beyond a critical value of A, the system enters the sym- 
metry breaking MQST phase, as predicted by the fixed point 
analysis earlier One of the novel aspects of the binary mix- 
ture is the existence of Zero-mode MQST accompanied by 
phase separation of the two species. Even with the initial con- 
ditions corresponding to both species abundance in the same 
well, the two components localize in the two different wells. 
In this case, transition to MQST is accompanied by phase sep- 
aration: although the two species overlap for some time, the 
< Za{t) > and the < Zf,(t) > have opposite signs. 



If the dynamics is initiated with 4>a {t = 0) = (pt {t = 0) = 
TT, both species oscillate in 7r-mode provided A < A° « 0.67, 
and initial population imbalance is small ( |Z;(0)| << 1) . 
Analogous to the Zero-mode, the dynamics in the tt mode is 
in general quasiperiodic. As seen in the figure, the motion is 
in phase with the slow mode and out of phase with the other. 
Comparison with the Zero and the 7r-mode oscillations show 



D. Symmetry Restoring and Phase Separation: 
Swapping-Mode 

As A increases further,the system exhibits "swapping- 
modes" where the two species swap places between the two 
wells but remain phase separated as shown in Fig. 4. As seen 
in the figure (at t — 0), the dynamics is initiated with positive 
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I (A=1.6. 1.8) (A=1.6. 1.8) 

FIG. 4: (color online) Transition from Josephson oscillations (top 
with A = 1.6 ) to MQST with phase separation (bottom, with A — 
1.8 ), obtained by varying A. The left and right plots show the time 
series and phase portraits respectively. 



population imbalance of both species. However, the result- 
ing dynamics corresponds to back and forth motion where the 
two species swap places between the two wells. In contrast 
to MQST, the swapping dynamics restores the symmetry of 
the tunneling solution in the double well. However, the two 
species remain mostly phase separated, avoiding each other 
by swapping. 

In other words, the swapping phase is characterized by < 
Zait) >=< Zb{t) = 0, but < Za{t)Zi,{t) < > . That is, 
at a given instant of time, the two species are more likely to 
be found in the separate wells. Thus in the swapping mode, 
the two species oscillate back and forth between the two wells 
and still manage to avoid each other The swapping is found 
to occur in the nonlinear Zero-mode as well as in the running 
mode . Furthermore, a transition from MQST to swapping 
phase can be achieved either by varying A ( Fig. 5) or by 
varying the initial conditions ( Fig. 6). 



E. Symmetry Breaking in 7r-Modes: Coexistence Phase 



For A < AJ^ w 0.67, and (t)a{t = 0) = (t)b{t = 0) = tt, 
and \Zi{t = 0)1 << l,both species execute small amplitude 
oscillations with < Zi{t) >— > and < (j)i{t) >~ tt, as 
shown in Fig. 7 Such modes are characterized by superpo- 
sition of sinusoidal modes with two competing frequencies 
and the resulting dynamics is in general quasiperiodic. As ex- 
pected from the fixed point analysis, the two species with both 
inter and intra-species repulsive interaction can self-trap in the 
same well. That is , we have MQST where the species coex- 
ist in the same potential well, inspite of repulsive interaction 
among them. 



O 50 lOO 150 200 




t (A = 2.3, 2.5) 

FIG. 5: (color online) Symmetry restoring transition by varying A 
where the upper panel with A — 2.3 shows MQST phase with phase 
separation while the lower panel with A = 2.5 shows phase separa- 
tion due to swapping mode. 




FIG. 6: (color online) Symmetry restoring transition obtained by 
changing initial conditions (Zi{t = 0)) slightly for fixed A = 2. 
Figure shows the time series as well as the phase portraits where the 
straight line and crosses show the corresponding fixed points. 



F. Swapping in 7r-modes 



As illustrated in figure [H within the 7r-mode phase, if the 
dynamics of the two species is initiated in separate wells, that 
is, Za{t = 0) and Zi,{t = 0) have opposite signs, the MQST 
phase can be destroyed when the initial population imbalance 
increases beyond a critical value. The tunneling solutions be- 
come symmetric as MQST is replaced by swapping modes. 
In this case the swapping can be viewed as the two species 
"chasing" each other 

It should be noted that the swapping dynamics in the Zero 
and the 7r-modes is very similar. However, swapping in the 
Zero-mode corresponds to two species avoiding each other 
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FIG. 7: (color online) Transition to MQST in tt modes. A = 0.6 
(top) and A — 0.8 (bottom) describe, respectively, the small ampli- 
tude TT-mode oscillations and MQST in vr-mode. 




V. EXPERIMENTAL REALIZATION 

The effects described in this paper should be realizable 
for condensate mixtures that already exist in the laboratory. 
One example in particular is a mixture of ^^Rb and ^^Rb 
atoms that has been created in several recent experiments at 
JILA iItI [12I1 . This system is relevant to the analysis in this 
paper because the scattering length, ass-ss. that characterizes 
the interaction between ^^Rb atoms is tunable by an exter- 
nal magnetic field via a feshbach resonance centered at ap- 
proximately 155 Gauss ifisll . Additionally, the interspecies 
scattering length, Oss-s?, is also tunable with two feshbach 
resonances (for a |2, — 2)85/|l, —1)87 collision) located at ap- 
proximately B = 267 Gauss and B = 356 Gauss. 
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FIG. 8: (color online) Transition from broken symmetry (MQST in tt 
modes) to symmetric configurations, obtained by changing the initial 
population imbalance. The three plots correspond to three different 
initial conditions. 



while swapping in the 7r-mode coiTesponds to one component 
chasing the other This is because, in the Zero-mode, species 
prefer residing in the separate wells while in the 7r-mode, they 
like to stay in the same well. This unique type of coherence 
between the two different species is one of the most fasci- 
nating aspect of the binary mixture dynamics in double well 
potential. 



In the most recent experiment UJ], a ^^Rb/^''Rb BEC 
mixture was produced by trapping a thermal-gas sample of 
the mixture and performing evaporative cooling on the ^^Rb 
which sympathetically cools the *^Rb. The cold gas mixture 
is then transferred to an optical trap that provides tight con- 
finement transverse to the trapping beam and loose confine- 
ment along the beam. If an additional pair of beams were 
applied along this direction as was done in the Albiez experi- 
ment list], it would create a setup to which the analysis in this 
paper would apply. 



VL SUMMARY 

Existence of a variety of BEC species with tunable inter and 
intra-species scattering lengths makes BEC mixtures one of 
the most attractive candidates for exploring novel phenomena 
involving quantum coherence and nonlinearity. Our analysis, 
based on the two-mode GP equation for the two interacting 
species of BEC in a double well trap unveils a variety of phe- 
nomena describing broken symmetry as well as subsequent 
restoration of symmetry, as we change the parameters or the 
initial conditions. Such coherence is found to exist over a 
broad range of parameters, establishing the robustness of the 
effects. 

To make direct comparison with experiments, we need to 
solve the coupled GP equations to obtain various parameters 
of the effective coupled pendulum system in terms of the mi- 
croscopic parameters of the system and work in this direction 
is in progress. Furthermore, by quantizing the Hamiltonian 
(coupled pendulum or the spin Hamiltonian), we hope to study 
quantum dynamics of number fluctuations that may code the 
emergence of new quantum phases in the system. 
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APPENDIX: TWO-MODE EQUATION PARAMETERS 

With (jx = g^N/h { X — a,b,ab ), the various coupling 
constants in the coupled equations can be shown to be given 

by, 

^lab = -gab J [(X-X-)' - ixUX)']dr 
Mab = -gab J [iX-X+f - ixlxtfjdr 

Aab = 2gab J {xU-xU-)dr 

Kb = [AE- fbAjb- faDab]/n 

Ca = (7++7- -27a)/2?i 

Cb = (7++7- -27"6)/2?i 

Dab = {A-Zab- A%b)/2h 



For each species, the localized spatial modes x± {x) are 
obtained by adding and subtracting the symmetric and the an- 
tisymmetric solutions of the time-independent coupled GPE 
equations. The AE is the difference in the chemical potential 
between the (symmetric) ground and the (anti-symmetric) first 
excited state of the coupled time-independent GPE equations. 



